Origin of the anomalous long lifetime of C 
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We report the microscopic origins of the anomalously suppressed beta decay of 14 C to 14 N using the 
ab initio no-core shell model (NCSM) with the Hamiltonian from chiral effective field theory (EFT) 
including three-nucleon force (3NF) terms. The 3NF induces unexpectedly large cancellations within 
the p-shell between contributions to beta decay, which reduce the traditionally large contributions 
from the NN interactions by an order of magnitude, leading to the long lifetime of 14 C. 

PACS numbers: 21.30.Fe, 21.60.Cs, 23.40.-s 



The measured lifetime of 14 C, 5730±30 years, is a valu- 
able chronometer for many practical applications ranging 
from archeology to physiology. It is anomalously long 
compared to lifetimes of other light nuclei undergoing the 
same decay process, allowed Gamow- Teller (GT) beta- 
decay, and it poses a major challenge to theory since tra- 
ditional realistic nucleon-nucleon (NN) interactions alone 
appear insufficient to produce the effect 1]. Since the 
transition operator, in leading approximation, depends 
on the nucleon spin and isospin but not the spatial coor- 
dinate, this decay provides a precision tool to inspect se- 
lected features of the initial and final states. To convinc- 
ingly explain this strongly inhibited transition, we need 
a microscopic description that introduces all physically- 
relevant 14-nucleon configurations in the initial and final 
states and a realistic Hamiltonian that governs the con- 
figuration mixing. 

We report the first no-core solutions of 14 C and 14 N 
using a Hamiltonian with firm ties to the underlying the- 
ory of the strong interaction, Quantum Chromodynam- 
ics (QCD), which allows us to isolate the key canceling 
contributions involved in this beta decay. We find that 
the three-nucleon force (3NF) of chiral perturbation the- 
ory (ChPT) plays a major role in producing a transition 
rate that is near zero, needed for the anomalous long 
lifetime. A chiral 3NF with coupling constants consis- 
tent with other works and within their natural range can 
provide the precise lifetime. This indicates that correc- 
tions to the lifetime that arise from increasing the basis 
space, from including additional many-body interactions 
and from corrections to the GT operator in ChPT [2[ 
may be absorbed into an allowed choice of the 3NF. 

Our work features two major advances over recent al- 
ternative explanations [H, H|: (1) we treat all nucleons on 
the same dynamical footing with the no-core shell model 
(NCSM) [|, and (2) we include the 3NF of ChPT || as 
a full 3-nucleon interaction. This follows previous work 
detailing the structure and electroweak properties of se- 
lected A=10-13 nuclei Q with the same chiral NN + 
3NF. We also establish a foundation for future work on 



the GT transitions to excited A=14 states [8|. 

ChPT provides a theoretical framework for inter- 
nucleon interactions based on the underlying symmetries 
of QCD. Beginning with pionic or the nucleon-pion sys- 
tem Q one works consistently with systems of increas- 
ing nucleon number [Iol - [l2l |. One makes use of the ex- 
plicit and spontaneous breaking of chiral symmetry to 
expand the strong interaction in terms of a generic small 
momentum. The ChPT expansion conveniently divides 
the interactions into perturbative and non-perturbative 
elements. The latter are represented by a finite set of 
constants at each order of perturbation theory that are 
not presently obtainable from QCD, but can be fixed by 
experiment. These constants should be of order unity 
("naturalness"). Once determined, the resulting Hamil- 
tonian predicts all other nuclear properties. We have 
demonstrated that this reductive program works well for 
light nuclei when the 3NF is included. 
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We adopt NN and 3NF potentials of ChPT 
and the no-core shell model (NCSM) to solve the many- 
body Schrodinger equation for A=14 nuclei while pre- 
serving all Hamiltonian symmetries 0, [TBI . flij . The non- 
perturbative coupling constants of the 3NF, not fixed by 
7r — N or NN data, are cd for the N — ir — NN contact 
term and ce for the NNN contact term. We previously 
fit cd and ce to A=3 binding energies and selected s-shell 
and p-shell nuclear properties and showed the resulting 
interactions worked well for A=10-13 nuclei Q- Here, 
we present results for (cd,Ce) = (—0.2, —0.205) and for 
(—2.0, —0.501). Both sets, labeled by their crj value be- 
low, are allowed by the naturalness criterium and fit the 
A=3 binding energies. The former also produces a pre- 
cise fit to the triton half life [17( . The latter produces the 
triton half life within 20% of experiment but is preferred 
by the 14 C half life as we show here. 

We select a basis space of harmonic oscillator (HO) 
single-particle states specified by the HO energy HQ 
and by N max , the limit on the number of HO quanta, 
^j(2nj + li), above the minimum for that nucleus, rii is 
the principal quantum number of the orbital of the i th 
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nucleoli and U is its orbital angular momentum. Thus 
Nmax — is the smallest basis space for each nucleus con- 
sistent with the Pauli principle. For each choice of N max 
and h£l we carry out a well-established finite-basis renor- 
malization 1 8] of the Hamiltonian to arrive at an Hermi- 
tian effective 3-body Hamiltonian [l6[ . All symmetries 
of the underlying Hamiltonian are preserved throughout. 
We adopt hVL — 14 MeV, which corresponds to the min- 
imum in the binding energy obtained for N max — 8. We 
have checked that results for 13 < Ml < 18 MeV do not 
qualitatively alter our conclusions. 

Basis space dimensions grow rapidly with N max and 
provide a major technical challenge [19|. The many-body 
matrix dimension in the M-scheme basis (total magnetic 
projection quantum number is fixed) for 14 C ( 14 N) is 
872,999,912 (1,090,393,922) at N max = 8 and the asso- 
ciated number of non- vanishing 3NF matrix elements on 
and below the diagonal ~ 2.9 x 10 13 (~ 3.9 x 10 13 ). 

We obtained our results on the Jaguar supercomputer 
at Oak Ridge National Lab [22| using up to 35,778 hex- 
core processors (214,668 cores) and up to 6 hours of 
elapsed time for each set of low-lying eigenvalues and 
eigenvectors. The number of non-vanishing matrix ele- 
ments exceeded the total memory available and required 
matrix clement recomputation "on-the-fly" for the itera- 
tive diagonalization process (Lanzos algorithm). 
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FIG. 1: (Color online) The excitation spectrum of 14 N ob- 
tained in the NCSM using chiral interactions as a function of 
basis space cutoff N max , the number that labels the columns. 
(J rlr ,T) values are in the legend. The center column displays 
the experimental spectrum for 14 N. Columns to the right are 
obtained with the chiral NN interaction alone while those 
on the left are obtained with the chiral NN + 3NF using 
cd = —0.2 (solid lines) and cd = —2.0 (dashed lines). 



interaction alone ("N3LO" columns) and the chiral NN 
+ 3NF ("N3LO+3NF" columns). 

The trends in excitation energies with increasing N max 
suggest reasonable convergence. We note that the level 
orderings are the same with and without the 3NF in the 
largest basis spaces {N max — 8). Overall, the spectra are 
more compressed with NN alone than with NN + 3NF, 
a pattern seen before @, EH. However, the 3NF appears 
to have over compensated for the compression in the NN 
spectra when compared with experiment. On the other 
hand the 3NF improves the binding energies. The exper- 
imental ( 14 C, 14 N) binding energies are (105.28, 104.66) 
MeV. Without the 3NF, we obtain (97.20, 96.36) MeV. 
When we include the 3NF the binding energies increase to 
(108.43, 108.47) MeV for c D = -0.2 and (108.30, 108.41) 
MeV for c D = -2.0 

We note that the binding energies for both 14 C and 
14 N are about 8 MeV or 8% larger at N max = 6 than 
they are at N max = 8. Similar calculations for lighter 
nuclei show that the convergence of the binding energy is 
not monotonic, hampering an extrapolation to the com- 
plete (infinite-dimensional) basis. On the other hand, 
the excitation spectra appear rather stable upon reach- 
ing N max = 8. We take results in the largest available 
basis space as estimates of the converged results. We also 
show in Fig. [T]that the spectra are relatively insensitive 
to a significant range of values for (cd, eg). 

Using Fermi's Golden rule for the transition rate, the 
half life T-^i% for 14 C is given by 



Tl/2 = 



27r 3 ?i 7 ln2 



f(Z,E ) m\c±G\ g A 2 \M GT \ 



(1) 



where Mqt is the reduced GT matrix element; f(Z,Eo) 
is the Fermi phase-space integral; Eq = 156 keV is the /3 
endpoint; Gy = 1.136 10~ n MeV~ 2 is the weak vector 
coupling constant; — 1-27 is the axial vector coupling 
constant; and m e is the electron mass. Mqt for the tran- 
sition from the initial 14 C (J*', T) = (0 + , 1) ground state 
to the 14 N (l+,0) ground state (#/) is defined by 
the spin-isopin operator a(k)T + (k) acting on all nucleons, 
k: 



(2) 



Since the initial state has total spin zero, Mqx is equal 
to the M-dependent GT matrix element MqL, 



a\(TT+\P) p a j3, 



(3) 



Fig. Q] compares our spectra for 1 N with experiment. 
The T = 1 isobaric analog states serve as good proxies 
for the 14 C spectra. We present the spectra as a function 
of the N max value (column labels) for both the chiral NN 



where (a|crT + |/3) is the one-body matrix element be- 
tween HO single-particle states a and /?, which is non- 
vanishing only when both single-particle states are in 
the same shell, and the one-body density matrix p a p = 

<*/I4o/j|*<). 
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In order to reproduce the measured half life of Ti/ 2 — 
5730 years, the GT matrix element must be anomalously 
small, I-Mq^I ~ 2 x 10~ 3 , in contrast with a conventional 
strong GT transition in a light nucleus with Mgt — 1. 
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FIG. 2: (Color online) Contributions to the 14 C beta decay 
matrix element as a function of the HO shell in which they are 
evaluated when the nuclear structure is described by the chiral 
interaction. Top panel displays the contributions with (two 
right bars of each triplet) and without (leftmost bar of each 
triplet) the 3NF at N max = 8. All contributions are summed 
within the shell to yield a total for that shell. The bottom 
panel displays the running sum of the GT contributions over 
the shells included in the sum. Note the order-of-magnitudc 
suppression of the p-shell contributions arising from the 3NF. 



In Fig. [21 we present our main results for Mqt- We 
decompose Mqt at N max = 8 into the contributions 
arising from each HO shell for two cases with the 3NF 
(cd = —0.2,-2.0) and one without. The largest effect 
occurs in the p-shell where, for both values of cd, the 
3NF reduces the contributions by an order-of-magnitude 
from the result with the NN interaction alone. The con- 
tributions of each of the 9 additional shells is enhanced by 
up to a factor of 2 by the 3NF. The cumulative contribu- 
tions to the GT matrix element is displayed in the lower 
panel of Fig. [2]where one sees clearly the net suppression 
to M G t due to the 3NF. 

Looking into the detailed changes within the p-shell, 
one finds that the 3NF introduces a systematic shift of 
strength away from 1-body density matrix terms involv- 
ing transitions between the 0p 3 / 2 and the 0pi/2 orbits 
to 1-body density matrix elements involving transitions 
within the Opi/2 orbits. The shift is about 30% of the 
magnitudes of these 1-body density matrix terms and 



supports a recurring theme of 3NF's in p-shell nuclei - 
they play a significant role in the spin-sensitive proper- 
ties of spin-orbit pairs. We note that our observed shifts 
within the p-shell due to the 3NF is similar to what Ref. 
[U accomplished with a density-dependent effective NN 
interaction obtained by modeling leading contributions 
from the chiral 3NF. However, our net contributions from 
other shells, which are absent in Ref. Q, overwhelm the 
net p-shell contribution as seen in Fig. [5] That is, while 
the s-shell and sd-shell contributions nearly cancel, all 
the shells above the sd-shell contribute about a factor 2 
greater than the now-suppressed p-shell contribution. 

To further understand the role of the 3NF, we can ex- 
amine the contributions to Mqt in the LS-scheme where 
the single-particle quantum numbers now involve the or- 
bital angular momentum projection mi and spin projec- 
tion m s replacing the total angular momentum projection 
nij. This is a convenient representation for this transi- 
tion since mi and m s must be the same for incoming and 
outgoing single-particle states. For the p-shell contribu- 
tions, the resulting decomposition to LS-scheme yields 
results shown in Table |U Note that there is nearly per- 
fect cancellation between the mi — and m; = ±1 terms 
once the 3NF is included. 

Given the overall effects on Mqt by inclusion of the 
3NF as seen through Fig. [5] and in Table HI one may fine 
tune cd to reproduce the value 2 x 10~ 3 consistent with 
the 14 C lifetime. This is analogous to the fine tuning of 
co in Ref. 17] to fit the 3 H lifetime. Our estimate that 
the desired suppression of Mqt occurs with crj — —2.0 
illustrates this point. Note that the spectra of 14 C and 
14 N are rather insensitive to this range of cd values. We 
show the resulting Mqt in the final row of Table HI 
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NN only 


NN + 3NF 
c D = -0.2 


NN + 3NF 
c D = -2.0 


(1>+*) 


0.015 


0.009 


0.009 




-0.176 


-0.296 


-0.280 


(o,4) 


0.307 


0.277 


0.283 


(o,-|) 


0.307 


0.277 


0.283 


(-i,+f) 


-0.176 


-0.296 


-0.280 


(-i,4) 


0.015 


0.009 


0.009 


Subtotal 


0.292 


-0.019 


0.024 


Total Sum 


0.275 


-0.063 


-0.013 



TABLE I: Decomposition of p-shell contributions to Mgt in 
the LS-scheme for the beta decay of 14 C without and with 
3NF. Lhe 3NF is included at two values of cd where cd — 
—0.2 is preferred by the 3 H lifetime and cd — —2.0 is preferred 
by the 14 C lifetime. The calculations are performed in the 
N m ax = 8 basis space with HQ, = 14 MeV. 



Next, we studied the sensitivity at (cd,ce) — 
(—2.0, —0.501) and at N max — 6 by successively setting 
each to zero while keeping the other fixed. The larger 
effect on Mqt appears with cd = 0. However, the re- 
sulting shifts in the magnitude of Mqt are approximately 
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proportional to the magnitude of the changes in cb and 
Ce which implies Mqt has about the same sensitivity to 
each. This sensitivity differs from that of Ref. Q where 
the Ce term was found to play a leading role in Mqt- 

Since both appear natural, we conclude that the dif- 
ferent values of cd from the beta decays of 3 H and 14 C 
accommodate the sum of higher-order effects: meson- 
exchange currents, additional ChPT interaction terms, 
and contributions from larger basis spaces. Additional 
3NF terms, e.g. those arising at the N 3 LO level of ChPT 
that do not impact the 3 H lifetime [2(|, and/or four nu- 
cleon interaction terms may impact 14 C and 14 N. 

In order to further test the ab initio wavefunctions and 
the effects of the 3NF, we present in Table HI] the con- 
tributions to selected electromagnetic properties of 14 N. 
We see that the 3NF has substantial influence on the 
magnetic-dipole transition in 14 N from the ground state 
(GS) to the (0+, 1) analog state of the 14 C GS though the 
mechanism for this suppression is quite different from the 
GT transition. In this Ml transition, we find significant 
cancellation between the nucleon spin and proton angu- 
lar momentum contributions. Without the 3NF, spin 
and angular momentum contributions add constructively. 
Electric properties such as the charge radius (RMS) and 
the electric quadrupole moment (Q) are reduced when we 
include 3NF due, in part, to the increased binding. Also, 
the GS magnetic moment (//) is modified by 10% when we 
include 3NF. We expect RMS and Q to increase with ba- 
sis space size while the GS energy, /i and the B(M1) will 
be less effected. The role of the 3NF on the electroweak 
properties of these nuclei is indeed multi-faceted. 



Observable 


Experiment 


NN only 


NN + 3NF 


NN + 3NF 
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c D = -0.2 


cd = -2.0 


RMS 


2.42(1) 


2.28 


2.25 


2.24 


Q 


1.93(8) 


1.87 


1.03 


1.19 




0.404 


0.379 


0.347 


0.347 


B(M1) 


0.047(2) 


1.002 


0.037 


0.098 



TABLE II: Properties of 14 N without and with 3NF in the 
Nmax = 8 basis space with ftfi = 14 MeV. The magnetic 
moments, which tend to converge rapidly, are obtained at 
N m ax = 6. The point proton root-mean-square radius (RMS) 
is quoted in fm. We corrected the measured charge radius 
(2.56(1) fm) for the finite proton charge contribution. The 
magnetic moment is in nuclear magnetons eh/2mc; and the 
quadrupole moment is in e 2 fm 4 (all for the GS). The B(M1) 
is the transition from the GS to the (0 + , 1) state (the isobaric 
analog of the 14 C GS). 

In conclusion, the chiral 3NF in ab initio nuclear 
physics produces a large amount of cancellation in the 
matrix element Mqt governing the beta decay of 14 C. 
This cancellation signals a major signature of 3NF ef- 
fects in the spin-isopin content of the Op orbitals in 14 C 
and 14 N. The 3NF, particularly through its longest range 
two-pion exchange component, modifies the Op spin-orbit 
pairs in the ground states of these nuclei to produce the 



suppression of Mqt consistent with the anomalous long 
lifetime of 14 C. 
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